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We study the nature of tunneling phase time for various quantum mechanical structures 
such as networks and rings having potential barriers in their arms. We find the generic 
1 presence of Hartman effect, with superluminal velocities as a consequence, in these sys- 

f**^ ' tems. In quantum networks it is possible to control the 'super arrival' time in one of the 

' arms by changing the parameters on another arm which is spatially separated from it. 

, This is yet another quantum nonlocal effect. Negative time delays (time advancement) 

and 'ultra Hartman effect' with negative saturation times have been observed in some 
parameter regimes. In presence and absence of Aharonov-Bohm (AB) flux quantum rings 
^ ' show Hartman effect. We obtain the analytical expression for the saturated phase time. 

In the opaque barrier regime this is independent of even the AB flux thereby generaliz- 
^ I ing the Hartman effect. We also briefly discuss the concept of "space collapse or space 

destroyer" by introducing a free space in between two barriers covering the ring. Further 
we show in presence of absorption the reflection phase time exhibits Hartman effect in 
contrast to the transmission phase time. 
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1. Introduction 

Quantum tunneling, where a particle has finite probability to penetrate a classically 
forbidden region is an important feature of wave mechanics. Invention of the tunnel 
diode^, the scanning tunneling microscope^ etc. made the quantum tunneling useful 
from a technological point of view. In 1932 MacColP pointed out that tunneling is 
not only characterized by a tunneling probability but also by a time the tunneling 
particle takes to traverse the barrier. There has been considerable interest on the 
question of time spent by a particle in a given region of space"'''^'^. The recent devel- 
opments of nanotechnology brought new urgency to study the tunneling time as it 
is directly related to the maximum attainable speed of nanoscale electronic devices. 
In addition, recent experimental results claiming superluminal tunneling speeds for 
photons^'^ call for detailed analysis of this problem. In a number of numerical^, 
experimental^' and analytical study of quantum tunneling processes, various def- 
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initions of tunneling times have been investigated. These different time scales are 
based on various different operational definitions and physical interpretations. Till 

date there is no clear consensus about the existence of a simple expression for this 
time as there is no hermitian operator associated with if*. Furthermore, a corpus- 
cular picture of tunneling is very hard to be realized due to the lack of a direct 
classical limit for the trajectories and velocities of the tunneling particle. 

Among the various time scales, 'dwell time^^^ which gives the duration of a 
particle's stay in the barrier region regardless of how it escapes can be calculated 
as the total probability of the particle inside the barrier divided by the incident 
probability current. The 'conditional sojourn time'^ gives the time of sojourn (dwell) 
in the spatial region of interest for some given conditions of scattering. It can be 
defined meaningfully by a 'clock' which is basically an extra degree of freedom 
that co-cvolves with the sojourning particle. Biittiker and Landauer proposed*^ 
that one should study 'timncling time' using the transmission coefficient through 
a static barrier of interest, supplemented by a small oscillatory perturbation. This 
is referred to as ' Biittiker-Landauer time^ in the literature. In another approach, 
the traversal time, is measured by the spin precession of the tunneling particle in 
a uniform infinitesimal magnetic field. This is called local 'Larmor time'^^'^^. A 
large number of researchers interpret the 'phase time'^'^'^'^^ as the temporal delay 
of a transmitted wave packet. This time is usually taken as the difFcrcncc between 
the time at which the peak of the transmitted packet leaves the barrier and the 
time at which the peak of the incident quasi-monochromatic wave packet arrives 
at the barrier. Within the stationary phase approximation, the phase time can be 
calculatcid from the energy derivative of the 'phase shift' in the transmitted or 
reflected amplitudes. It would be worth mentioning that for a symmetric barrier, 
the reflection and transmission phase times are equal. For tunneling through a ID 
static barrier due to the time-reversal symmetry, one can show from the unitarity of 
the scattering matrix that the phases of the reflection and transmission amplitudes 
differ by a constant, quantitatively by 7r/2. However this does not remain true for 
asymmetric/complex barrier. Biittiker-Landauer*^ raised objection that the peak 
is not a reliable characteristic of packets as the wave packet may undergo strong 
distortion or deformation after tunneling through the barrier. Moreover, there is no 
causal relationship between the peak of the transmitted wave packet and the peak of 
the incident packet. This is due to the fact that the peak of the transmitted packet 
can leave the scattering region before the peak of the incident packet has arrived. In 
contrast to 'dwell time' which can be defined locally, the 'phase time' is essentially 
asymptotic in character*^. The 'phase time' statistics is intimately connected with 
dynamic admittance of micro-structures''^''. The 'phase time' is also directly related 
to the density of states^*' The universality of 'phase time' distributions in 
random and chaotic systems has already been established earlier^^. In the case of 
'not too opaque' barriers, the tunneling time evaluated either as a simple 'phase 
time'^ or calculated through the analysis of the wave packet behaviour^'^ becomes 
independent of the barrier width. This phenomenon is termed as the Hartman 
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effect 15'23.24 xhis implies that for sufficiently long barriers the effective velocity 
of the particle can become arbitrarily large, even larger than the light speed in the 
vacuum (supcrluminal effect). Though this interpretation is a little far fetched for 
non-relativistic Schrodinger equation as velocity of light plays no role in it, this 
effect has been established even in relativistic quantum mechanics. 

Tunneling time associated with electron being very small (-^ femto second) the 
experiments to verify the 'Hartman effect' (generalized Hartman effect) are usu- 
ally done on optical waveguide where the corresponding time is of the order of 
pico second. However the formal analogy between the Schrodinger equation and the 
Helmholtz cqiiation for electromagnetic wave enables one to correlate the results 
for optical experiments to that for electrons. Photonic experiments show that elec- 
tromagnetic pulses travel with group velocities in excess of the speed of light in 
vacuum as they tunnel through a constriction in a waveguide^. Experiments with 
photonic band-gap structures clearly demonstrate that 'tunneling photons' indeed 
travel with supcrluminal group velocities^. Their measured tunneling time is prac- 
tically obtained by comparing the two peaks of the incident and transmitted wave 
packets. Thus all these experiments directly or indirectly confirmed the occurrence 
of Hartman effect without violating 'Einstein causality' i.e., the signal velocity or 
the information transfer velocity is always bounded by the velocity of light. It should 
also be noted that in the photonic tunneling time experiments by Nimtz et. al, based 
on frustrated total internal reflection, the velocity of the half-width of the pulse (not 
the peak of the wave packet) is monitored. The velocity of the half-width is found 
to be supcrluminal (according to theory and experiment). Relation of this result to 
the causality principle is discussed in the references'^ The 'Hartman effect' has 
been extensively studied both for nonrelativistic (Schrodinger equation) and rela- 
tivistic (Dirac equation)^'^'^ cases. Recently Winful showed that the saturation 
of phase time is a direct consequence of saturation of integrated probability density 
under the barrier (cquivalently in the electromagnetic waves saturation of stored 
energy) . The Hartman effect has been found in one dimensional barrier tunneling'^ 
as well as in tunneling through mesoscopic rings in presence of Aharonov-Bohm 
(AB) flux'*. In the present work we extend the study of phase times for branched 
quantum networks and rings. 

The main results of this paper are as follows: 

• We have found Hartman effect in barrier tunneling regime for quantum 
networks and rings. 

• In quantum networks non-locality and time-advancement (negative time 
delay) are found. 

• In presence of AB-flux in the rings saturated phase time becomes indepen- 
dent of the flux. Thus the Hartman effect is obtained even in presence of 
AB-flux. 

• We obtain space collapse in quantum ring i.e. in presence of inter-barrier 
free space the phase time becomes independent of the length of this free 
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space in the off resonance regimes. 
• Even in presence of absorption in the barrier, the reflection phase time 
shows Hartman eff'ect though the transmission phase time grows with the 

length of the barrier. 

The paper is organized as follows. In section 2 we present phase times^^ for 
branched networks of quantum wires which can readily be realized in optical wave 
propagation experiments. This geometry allows us to chock other nonlocality effect 
such as tuning of the saturation value of 'phase time' and consequently the super- 
luminal speed in one branch by changing barrier strength or length in any other 
branch, spatially separated from the former. In section 3 we study the Hartman 
effect on a quantum ring geometry i. e. beyond one dimension and in the presence of 
AB-flux^*''^°. Our results confirm 'Hartman effect' in quantum ring even in presence 
of AB-flux. Then we study the phase time in quantum ring in absence of AB-flux. In 
tunneling through the ring having two potential barriers with an intermediate free 
space where a quantum particle can propagate as a traveling wave, the saturated 
phase time becomes independent of the intermediate free length (in the large length 
limit of the barriers) in the off resonant cases. This result can be interpreted as a 
"space collapse or space destroyer" as if the intermediate free space does not ex- 
ist. We also show that in the presence of absorption in barrier regime the Hartman 
effect survives in the reflection mode as opposed to the transmission mode. Then in 
section 4 we summarize our results and conclude. 

2. Hartman effect and non-locality in quantum networks 

Hartman effect is itself one of the manifestations of quantum non- locality ^. Here we 
study the effect for various quantum mechanical networks having potential barriers 
in their arms. In such systems it is possible to control the 'super arrival' time in one 
of the arms by changing parameters on another, spatially separated from it. This is 
yet another quantum nonlocal effect. Negative time delays (time advancement) and 
'ultra Hartman effect' with negative saturation times have been observed in some 
parameter regimes. 

As a model system we choose a network of wires. The length of these wires are 
so narrow that only the motion along the length of the wires is of interest (a single 
channel case). The motion in the perpendicular direction is frozen in the lowest 
transverse sub-band. In a three-port Y-branch circuit (Fig. 1) the two side branches 
of quantum wire Si and are connected to a 'base' arm at the junction J. In 
general one can have A''(> 2) such side branches connected to the 'base' wire. 

We study the scattering problem across a network geometry as presented 
schematically in Fig. 1. Such geometries arc important from the point of view of 
basic science due to their properties of tunneling and interference''^''^^ as well as 
of technological applications such as wiring in nano-structures. In particular, the 
Y-junction carbon nanotubes are in extensive studies and they show various inter- 
esting properties like asymmetric current voltage characteristics"^^. In our system of 
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Fig. 1. Schematic diagram of a Y-junction or three-way spUtter. 

interest there are finite quantum mechanical potential barriers of strength Vi and 
V2 and length Ibi and /&2 in the side branches 6*1 and ^2 respectively. We focus on a 
situation wherein the incident electrons have an energy E less than Vn, n = 1,2. The 
impinging electrons in this sub-barrier regime travel as an evanescent mode/ wave 
and the transmission involve contributions from quantum tunneling and multiple 
reflections between each pair of barriers and the junction point. Here we are inter- 
ested in a single channel case where the Fermi energy lie in the lowest sub-band. To 
excite the evanescent modes in the side branches one has to produce constrictions 
by making the length of the regions of wires containing barriers much thinner than 
that of the other parts of the wire. The electrons occupying the lowest sub-band in 
the connecting wire on entering the constrictions experience a potential barrier (due 
to higher quantum zero point energy) and propagate as an evanescent mode^^'^^. 
In this work an analysis of the phase time or the group delay time in such a system 
is carried out. 



2.1. Theoretical treatment 

We approach this scattering problem using the quantum wave guide theory ^7,38^ 
the stationary case the incoming particles are represented by a plane wave 6**°^ of 
unit amplitude. The effective mass of the propagating particle is m and the energy 
is E = h^k"^ /2ra where k is the wave vector corresponding to the free particle. The 
wave functions, in different regions of the system considered in Fig. 1 can be written 
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as, 



Ipin (Xo) 

V'(2),(a;2) 
V'(i)7,(a;i) 
V'(2)„(a;2) 

V'(2),„(a;2) 



+ Re 



-ikxn 



(in5o) , 



^^gifexi ^^^g-ifexi (region I in Si), 
^2 ^= + e"' ^= (region I in ^s) , 
Ci e-'^i (^^-^^^ + Di e""' (^^-^^) (region II in5i) 
Ca 6"'== + £>2 e-^^ (region II in^a) 

■* (region III inSi) , 



to e 



i k {X2—W2 — Ib2) 



(region III inSa) , 



with K„ = ^2m(K - £)/fi2 being the imaginary wave vector in presence of a 
rectangular barrier of strength y„, with n = 1, 2. is the spatial coordinate for the 
'base' wire, whereas xi and X2 are the spatial coordinates for the and 5a side 
branches respectively. All these coordinates are measured from the jimction J. wi 
and W2 are the distances between the junction J and the respective barriers in arm 
Si and ^a- 

To solve the problem, we use Griffith's boundary conditions^® 



tpin{X0 = J)= V'(l)/(a;i = J) = tp{2)i{xi = J), 



and 



dtpin{Xo) 



dxo 



3X2 



(1) 



(2) 



J dxi 

at the junction J. All the derivatives are taken either outward or inward from the 
junction^^. In each side branch, at the starting and end points of the barrier, the 
boundary conditions can be written as 

'^(n)i{Xn=Wn)='>l^{n)ii{Xn=Wn), (3) 
{Xn =Wn + lhn)= ^(n)iii {Xn = Wn + Ihn) , (4) 



dx„ 



dXr. 



(u)„+!b„) ' 



(5) 
(6) 



(ii>„+i6„) dXn 

where n = 1,2. From the above mentioned boundary conditions one can obtain 
the complex transmission amplitudes ti and ia on the side branches Si and 5*3 
respectively. 



2.2. Results and Discussions 

Following the method introduced by Wigner^^, we can calculate the 'phase time' 
(phase time for transmission) from the energy derivative of the phase of the trans- 
mission amplitude tn^'^^ as 



dArg[t, 
dE 



(7) 
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Fig. 2. For a 3-way splitter with a barrier in Si arm, the 'phase times' ri and T2 are plotted as 
a function of barrier length Hbi' in (a) and (6) respectively. The solid, dotted, dashed curves are 
for Vi = 5, 4, 3 respectively. Other system parameters are E = l,wi = 3. 



In what follows, we set h = 1 and 2m = 1. We now proceed to analyze the 
behavior of t„ as a function of various physical parameters for different network 
topologies. We measure time at the far end of each barrier in the branched arms 
containing barriers and in the case of arms in absence of any barrier we measure 
the phase time at the junction points. We express all the physical quantities in 
dimensionless units i.e. all the barrier strengths Vn in units of incident energy 
E {Vn = Vn/E), all the barrier lengths Ibn in units of inverse wave vector 
[Ibn = klbn), where k = \fE and all the extrapolated phase time r„ in units of 
inverse of incident energy E (t„ = ETn). 

First we take up a system similar to the Y-junction shown in Fig. 1 in presence 
of a barrier Vi of length Ibi in arm but in absence of any barrier in arm 82- For 
a tunneling particle having energy E <Vi we find out the phase time n in arm Si 
as well as T2 in arm ^2 as a function of barrier length Ibi (Fig. 2). From Fig. 2(a) 
it is clear that Ti evolves with Ibi and eventually saturates to Ti^ for large Ibi to 
show the Hartman effect. Fig. 2(b) shows the phase time T2 in arm 52 which does 
not contain any barrier. This also evolves and saturates with Ibi, the length of the 
barrier in the other arm 5*1 . This delay is due to the contribution from paths which 
undergo multiple reflection in the first branch before entering the second branch 
via junction point J. In absence of a barrier in the n-th arm the phase time t„ 
measured close to the junction J should go to zero i.e. t„ ^ in the absence of 
multiple scatterings in the first arm. Note that tis and T2s change with energies of 
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Fig. 3. Here for a 3- way splitter witii one barrier in each branched arm Si and 5*2, the 'phase 
time' Ti is plotted as a function of barrier length '/6i' keeping lb2{= 1) and Vi(= 5) fixed and for 
different values of parameter V2- The dashed, solid and dotted curves are for V2 = 2.5,5.0 and 
12.5 respectively. Other system parameters are E = l,wi = W2 = S. In the inset ris is plotted as 
a function of V2 for the same system parameters. 



the incident particle (Fig. 2). From Fig. 2 it can be easily seen that T2s is always 
smaller than for any particular Vi i.e. the saturation time in the arm having no 
barrier is smaller. The phase time in both the arms show non-monotonic behavior 
as a function of Vi. As we decrease the strength of the barrier Vi the value of ti 
(T2) decreases in the whole range of lengths of the barrier and also the saturated 
value of Tis (t2s) decreases until Vi reaches 1.6 and with further decrease in Vi the 
values of Ti {T2) as well as Tis (t2s) starts increasing. 

As the second case we take up another Y-junction which contain potential bar- 
riers in both of its side branches as shown in Fig. 1. We fix the values of Vi{= 5) 
and vary Ibi for each values of V2 to study the Z6i-dependence of ri (Fig. 3). From 
Fig. 3 we see that ti decreases with increase in Ibi and saturate to a value tis at 
each value of V2 thereby showing 'Hartman effect' for arm '5*1'. But now, we can 
tune the saturation phase time i.e. speed of the peak of the wave packet in one arm 
Si non- locally by tuning the strength of the barrier potential V2 sitting on the other 
arm 6*2! Thus 'quantum nonlocality' enables us to control the 'super arrival' time in 
one of the arms (Si) by changing a parameter (V2) on the other, spatially separated 
from it. In the inset of Fig. 3 we plot rig as a function of V2. It clearly shows that 
when the barrier strengths Vi and V2 arc very close the 'phase time' reaches its 
minimum value. In all other cases i.e., whenever Vi ^ V2, the value of Tig is larger. 
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Fig. 4. Here for a 3- way splitter with an identical barrier in each branched arm Si and 52, 
the 'phase times' arc plotted as a function of barrier length 'lbi'{= Z62) for fixed Vi = V2(= 5). 
The solid and dashed curves represent ti (=T2) and Tr respectively. Other system parameters are 
E = l,wi = W2 = 2.5. In the inset Tia(= T2s) is plotted as a function of V{= Vi = V2) for the 
same system parameters. 



With the above system but for two identicaUy branched arms (i.e., Vi = V2, 
wi = W2 and Ibi = lh2) we study the behaviour of different phase times as a func- 
tion of Ibi. For this case we see from Fig. 4 that n and T2 are the same, as expected. 
The phase time in reflection mode, measured near the junction on base arm, evolves 
differently from these two phase times. The reflection phase time getting saturated 
with increasing barrier's length confirms Hartman effect. The value of saturated 
reflection phase time is also different from that of the transmission mode. We also 
study the saturated transmission phase times by changing the strength of the bar- 
riers. From the inset of Fig. 4, we see that for strength very close to the value of 
the incident energy the saturated phase time is quite large and the saturated phase 
time decreases with the increasing barrier strength. 

Wc shall now show another interesting result related to the Hartman effect. For 
this we keep V2(= 5) unaltered and reduce lh2 - For very small lb2{= 0.5) we see from 
Fig. 5 that t\ is negative for almost the whole range of Z6i-values showing 'time- 
advancement' and eventually after a sharp decrease saturates to a negative value 
of Tis = —4.514 implying 'Hartman effect' with advanced time or 'ultra Hartman 
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Fig. 5. Here for a .3-way splitter witli one barrier in eacli side branch Si and S2, the 'phase time' 
Ti is plotted as a function of Hbi' for a very small lb2{= 0.5). Other system parameters are E = I, 
wi = W2 = 2.5, V2 = 5 and Vi = 15. In the inset, the solid and dashed curves represent T2 and 
|t2p respectively as a function of Ibi. For better visibility we have plotted phase times scaled down 
by a factor of 20. 

effect'. In the inset we plot the eorresponding T2 and |i2p as a function of Ibi. The 
tunnehng phase time for an one dimensional barrier with the same strength V2 and 
width Z62, as used in the side branch ^2 of the network, is different from T2, obtained 
for the whole range of values of Ibi. This is due to the presence of the other arm 
in the network problem. In the cases discussed so far T2 vary more sharply in small 
Ibi regime. Furthermore, the inset in Fig. 5 shows a dip in T2 at parameter regimes 
where |f2p has a minimum. 

For a wave packet with large spread in real space it is possible that the leading 
edge of the wave packet reaches the barrier much earlier than the peak of the packet. 
This leading edge in turn can tunnel through the barrier to produce a peak in the 
other end much before the peak of the incident wave packet reaches the barrier 
region. This sometimes is referred to as pulse reshaping effect. This, in general, 
causes 'time advancement'^. This negative delay does not violate causality, however, 
this delay time is bounded from below. In ID barrier such a situation does not arise. 
In presence of square wells in one dimensional systems negative time delays have 
been observed. This effect is termed as 'ultra Hartman effect' [ see for dctails^*^ ]. 

As the next case we set the length of the barrier in arm 5*1 of the above mentioned 
system at a large value (say 100) where all the phase times get saturated. Now we 
shift the position of the barrier in arm 5*2 away from the junction and study its 
effect on the saturated transmission and reflection phase times. In Fig. 6 we have 
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Fig. 6. For a 3-way splitter with one barrier in each side branch Si and S2, ris/SO (thick soUd), 
T2s/50 (thick dotted) and rrs/50 (dashed) are shown as a function of '102'- |t2p (sohd) and |rp 
(double dot-dashed) are shifted upwards along the the positive y-direction by 3. The different 
system parameters are E = 1, Vi = 15, V2 = 5, Ibi = 100.0, lb2 = 0.5 and wi = 2.5. 



plotted all these three quantities ti,, and r^g as a function W2. Earlier we had 
shown in Fig. 3 that by changing a nonlocal parameter V2 one can tune tis whereas 
Fig. 6 shows that change of another parameter W2 can tune Tig nonlocally. Note that 
T2s and Trs also depend on 11)2- From Fig. 6 we see that 1^2^ shows resonances 
(anti-resonances) as a function of W2. The saturated delay time tis exhibit sharp 
variation in phase time around the resonance. Moreover, it takes negative as well 
as positive values. Around the resonance sharp peaks exist in both positive and 
negative side. In contrast to this observation of tis the other two phase times show 
only positive peaks. As |t2p and |rp have finite non-zero values, these variations 
in phase times should, in principle, be observable in experiments. In Fig. 7 we have 
plotted Tis as a function of W2 for two different values of the length Z62 of the 
barrier in arm 5*2. Note that as we increase the length Z62, the frequency of getting 
negative saturation values (ti^) reduces (see the solid curve in Fig. 7) and increasing 
the length Z62 further, the negative saturation goes away. This is in agreement with 
the discussions in previous paragraph. The saturation to a negative value can be 
attributed to the nonidcntical barriers in two branched arms, the barrier length of 
one being very small compared to the other. When Ibi = lb2 the time- advancement 
goes away. 

So far our discussion is based on wave guide theory of transport. We can as well 
treat the junction as a scatterer described by the following one parameter scattering 
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Fig. 7. Here for a 3- way splitter with one barrier in each side branch Si and S2, the 'saturated 

phase time' ri^ is plotted as a function of 'w2'. The dashed and solid curves are for lh2 = 0.5 and 
2.0 respectively. Other system parameters are E = I, Vi = 15, V2 = 5, Ibi = 100.0 and wi = 2.5. 



matrix*^ 

Sj=\ a b \ (8) 

V ha) 

where a = ^ — 2e — l) and 6=5 (Vl — 2e + l) . e is a coupling parameter with 
values < e < 0.5. When e — !■ the side branches are decoupled from the base arm 
while for e — > 0.5 they arc strongly coupled. This S'-matrix satisfies the conservation 
of current For our network system with two side branches, the equations at 
the junction J can be written as 



(9) 



We have studied the effect of the coupling parameter e on the Hartman effect. 
From Fig. 8, we see that though the nature of ti and T2 as a function of Ihi does 
not depend on the coupling parameter e but their saturation value is strongly de- 
pendent on e. For weak coupling, the value of saturated phase time is quite low 
and it increases with increasing strength of the coupling between base arm and side 
branches. For e = 4/9 we get back the results obtained earlier using the wave-guide 
theory (compare the solid curves in Fig. 8 with the dashed curves in Fig. 2). 
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Fig. 8. For a 3-way splitter with a barrier in Si arm, the 'phase times' ri and T2 are plotted as a 
function of barrier length 'Ibi' in (a) and (6) respectively. The dotted, dashed and solid curves are 
for the junction coupling parameter e = 1/9,1/3,4/9 respectively. Other system parameters are 
Vi = 3, E = l,wi =3. We see that both and T2s increase as we increase the coupling strength. 



3. Hartman effect in Quantum ring 

In this section wc study the scattering problem across a quantum ring connected to 
one ideal semi infinite lead (as shown schematically in Fig. 9). Such ring geometry 
systems are extensively investigated in mesoscopic physics in analyzing normal state 
Aharanov-Bohm effect which has been observed experimentally^^'^^. A magnetic 
field is applied perpendicular to the plane of the ring. Due to this a magnetic flux 
4> as shown in Fig. 9, enclosed by the ring, there is a finite quantum mechanical 
potential of strength V inside the ring while that in the connecting lead is set to 
be zero. We focus on a situation wherein the incident electrons have an energy 
E less than V. The impinging electrons in this sub-barrier regime travels as an 
evanescent mode throughout the circumference of the ring and the reflection or the 
conductance involve contributions from both the Aharanov-Bohm effect as well as 
quantum tunneling. Here we are interested in a single channel case where the Fermi 
energy lies in the lowest sub-band. To excite the evanescent modes in the ring we 
have to make the width of the ring much less than that of the connecting lead. The 
electrons occupying the lowest sub-band in the lead on entering the ring experience 
a higher barrier (due to higher quantum zero point energy) and propagate in the 
ring as evanescent mode. The transmission or conductance across such systems has 
been studied in detail^^'^^. An analysis of the phase time for such a ring system is 
carried out in the following subsections. 

Now we discuss the 'reflection phase time' for the above quantum ring in pres- 
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Fig. 9. Schematic diagram of a circular ring connected to semi infinite lead. 



ence of AB-flux as shown in Fig 9(a). It is well known that in one dimensional scat- 
tering/tunneling problem, reflection involves prompt part as well as the multiple 
scattering arising from the edges of the scattering center (say, for the square bar- 
rier). However, transmission across the scattering region docs not have the prompt 
part but has only contributions from multiple scattering. We would like to empha- 
size here the fact that the unitarity of the scattering matrix forces transmission 
and reflection phase times to coincide for a one dimensional tunneling problem (to 
be identical in magnitude), even though reflection has a prompt part as mentioned 
above. Hence, the information we get does not depend on whether we study the 
phase time in the reflection or in the transmission mode. Thus in the present sec- 
tion we have chosen a simple and generalized geometry where we can analytically 
study the phase time in a reflection mode and in presence of Aharonov-Bohm flux. 
We show that this phase time in the opaque barrier regime becomes independent of 
the length of the circumference of the ring and the magnitude of the AB-flux. We 
have also studied this effect by including an additional potential well between two 
barriers in the circular ring (Fig. 9 (b)). Interestingly, the saturated reflection phase 
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time becomes independent of the length of the potential well (in the large length 
limit) for energy away from resonances. Inside the potential well the electron trav- 
els with a real velocity. Increasing or decreasing the free path (length of the well) 
does not alter the saturated reflection phase time through the system. It seems as if 
for the electronic wave the free space collapses. This result is regarded as a "space 
collapse" or "space destroyer" 

3.1. Theoretical Treatment 

We use the same quantum wave guide theory^^'^^ as discussed in the earlier section. 
To set the equations in most general context we consider Fig. 9(b) in presence of 



flux (j). The wave functions in different regions are 

^o(xo) = e*''^" + re-**^^" ( in region I ) (10) 

ipi{xi) = Aie^ '''''' +316"' (in region II) (11) 

V'»(a;^) = Ce^'^^™ +L»e-*'=^™ ( in region III) (12) 

?A2(a;2) = ^2 6"'"""'" +^2 6'"=^"'" (in region IV) (13) 



with k being the wave-vector of electrons in the lead and in the intermedi- 
ate free space between two barriers inside the ring, m = ^/2m{Vi — E)/h'^ and 
K2 = a/2to(V2 — E)IW- are the imaginary wave-vectors respectively for tunneling 
electrons in the barriers of strength V\ and Vi inside the ring. The origin of the 
co-ordinates of xq and x\ is assumed to be at Ji and that for and xi at P\ and 
Pi respectively. At P\^x\= Ibi, at P2, x^ = w and at Ji, X2 = lb2, where Ibi and 
/&2 arc the length of the two barriers and w is the length of the well inside the ring. 
Total circumference of the ring is L = Ibi + lb2 + w. 

In presence of the AB-flux, following the same method described above, the 
boundary conditions for the current system (shown in Fig. 9(b)) are 

1 -Fr- Ai - Biexp(-iQ!i) = 0, (14) 
A2 exp{-K2lb2) exp(ia2) + B2 exp(K;2Z62) 

-l-r = 0, (15) 

ik (1 — r) + Ki (Ai — Bi cxp(— i ai)) — K2 A2 
exp(-K2 lb2) cxp(« ^2) - K2 B2 cxp(fi;2 lb2) = , (16) 
Ai exp(— Ki Ibi) exp[iai] 4- Bi cxp(Ki Ibi) 

-C-£)exp(-ia^) =0, (17) 
Ki Ai exp(— Ki lb\) exp(i ai) — /ti Bi exp(K;i Ibi) 

+ ikC - ikD ex.p{-iayj) = 0, (18) 
C exp{ikw) exp(iai„) -|- D exp(— ifcw) 

- ^2 - ^2 exp(-ia2) = , (19) 
ikC exp{ikw) — ikD exp{—ikw) exp(— ia^) 

-hK2 ^2 - K2-B2 exp(-ia2) = 0, (20) 
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Fig. 10. In absence of magnetic flux (i. e. (j> =0 ), for a ring with a barrier of strength V throughout 
its circumference, the reflection phase time Tr is plotted as a function of ring's circumference L. 
The solid, and dotted curves are for y = 5, 3 respectively. Incident energy is set to be S = 1. In 
the inset the saturated value of phase time Tr s is plotted as a function of the barrier's strength 
for same E. 



where iai, i (ai + a^) are the phases picked up respectively at Pi and P2 by the 
electron traveling clockwise from Ji and i (ai + a^, +a2) is the phase picked up by 
the same electron at Ji after traversing once along the ring. The total phase around 
the ring becomes ai + + 0:2 = '2 (t> / (t>o ■ 

To obtain an analytical expression for the reflection amplitude for a ring system 
as shown in Fig. 9(a) we solve Eqs. (14)- (20) using C = 0, D = 0, A2 = Ai, 
B2 = Bi, K2 = Ki in the wave functions (12) and (13) and a;i = L at Ji after a 
complete traversal along the circumference. The reflection amplitude is 

— Ki (2 cos(q:) - exp(A;L)) + i | exp(A:L) ^^^^ 
Ki (2 cos(a) — exp{kL)) + 1 1 exp{kL) ' 

where a = 2 7r0/0o- After knowing r, the 'reflection phase time' Tr can be calculated 
from the energy derivative of its phase^'^^ as 

dArg[r] 

= ^^. (22) 



I 



I 
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3.2. Results and Discussions 

We now proceed to analyze the behavior of Tr as a function of various physical 
parameters for these ring systems. In the similar fashion, described above, here also 

we express all the physical quantities in dimcnsionlcss units. Thus the reflection 
phase time Tr is expressed in units of inverse of incident energy E (rr = Err). After 
straight forward algebra in the large length (L) limit and in absence of magnetic 
flux, wc obtain an analytical expression for the saturated reflection phase time (using 
Eq. (21) in Eq. (22)), which is given by. 



with K being the imaginary wave vector of the electron inside the barrier of strength 
V and length lb (= L). 

First we take up a ring system with a single barrier along the circumference of 
the ring as shown in Fig. 9(a). For a tunneling particle having energy less than the 
barrier strength we find out the reflection phase time r,. as a function of barrier 
length L which in turn is the circumference of the ring. We see (Fig. 10) that in 
absence of magnetic flux, Tr evolves as a function of L and asymptotically saturates 
to a value Trs which is independent of L thus confirming the 'Hartman effect'. From 
Fig. 10 it is clear that the saturation value increases with the decreasing barrier- 
strength. In the inset of Fig. 10, we have plotted « as a function of barrier-strength. 
From this we can see that for electrons with incident energy close to the barrier- 
strength the value of t,. is quite large. 

To see the effect of magnetic fiux on 'Hartman effect', we consider the same 
system but in presence of AB-flux. We study the reflection phase time as a function 
of AB-flux for different values of the length L of the ring's circumference. We have 
chosen the lengths such that in the absence of the 'AB-flux', for a given system (i.e., 
for known E and V) the reflection phase time Tr gets saturated for these lengths. 
From Fig. 11 we see that as a function of (p shows AB-oscillations with an average 
value which is the saturation value Tr s for the same system in absence of AB-flux. 
Further we observe that (Fig. 11) is flux periodic with periodicity (/>o- This is con- 
sistent with the fact that all the physical properties in presence of AB-flux across 
the ring must be periodic function of the flux with a period ^o^'^'^^'^^- However, we 
see that as we increase L the magnitude of AB-oscillation in decreases. Conse- 
quently in the large length limit the visibility of oscillations vanishes. This clearly 
establishes 'Hartman effect' even in presence of AB-flux. The constant value of Tr 
thus obtained in the presence of flux is identical to Trs (0.294115) in the absence of 
flux (see Fig. 11) in the large length regime and its magnitude is given by Eq. (23). 
This result clearly indicates that the reflection phase time in the presence of opaque 
barrier becomes not only independent of the length of the circumference but also 
is independent of the AB-flux thereby observing the 'Hartman effect' even in the 
presence of AB-flux. Since the magnetic field tunes the boundary conditions and in 




(23) 
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Fig. 11. For a ring with a barrier of strength V lies throughout its circumference, the saturated 
phase time Tr s is plotted as a function of magnetic flux 4>. The soUd, dashed and dotted curves 
are for L = 6, 7, 9 respectively. Other system parameters are V = 5, E = 1. 

the largo length limit the evanescent states become insensitive to the boundaries 
(as they decay out before reaching the boundary) the Hartman effect follows. 

Now we consider the ring system with the ring having two successive barriers 
separated by an intermediate free space (we shall call it 'qiiantum well' in what 
follows) as shown in Fig. 9(b). In the absence of magnetic flux, we see the effect of 
'quantum well' on the reflection phase time Tr- In Fig. 12 Tr is plotted as a function 
of one of the barrier length (say Ibi) while other barrier length is kept fixed (/&2 = 5) 
and for few different values of length of the well. Here, the fixed value of lb2 is chosen 
in such a way that in absence of the well region the reflection phase time reaches 
saturation at this length. From Fig. 12 we see that for all parameter values of well 
length, the saturation value of reflection phase time Trs is same and it is equal to 
what we obtained in absence of the well in the ring system. Thus the saturated 
phase time becomes independent of the length of the well (in the long length limit) 
for the energy away from the resonances. This is as if the effective velocity of the 
electron in the well becomes infinite or equivalontly length of the well does not count 
(space collapse or space destroyer) while traversing the ring. 

Next we consider a similar system as that shown in Fig. 9(b). Here we see 
the effect of resonances present in the ring system with the well, on the saturated 
reflection phase time Trs- For the system described above with two identical barriers 
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Fig. 12. In absence of magnetic flux (i.e. (j> =0 ), for a ring with two barriers of strength Vi and 
V2 separated by an intermediate well region, the reflection phase time r,. is plotted as a function 
of ring's circumference L for different length w of the well. The dotted, solid and dashed curves 
are for «; = 0, 5, 10 respectively. Other system parameters are lb2 = 5, Vi = V2 = 2 and E = I. 



of strength Vi = V2 = 2 and length Ibi = lb2 = 5 separated by an intermediate free 
space of length w we study the phase time as a function of w in absence of the 
magnetic flux {(p = 0). In Fig. 13(a) we have plotted r^s, for the electrons with an 
incident energy E = 1.5. It should be noted that as we increase the length of the 
well the incident energy E coincides with resonances (or resonant states) in the well 
(which arise due to constructive interference due to multiple scatterings inside the 
well). For these values of w we observe sharp rise in the saturated delay time and 
its magnitude depends on the length of the well. We see that the resonances which 
have Lorentzian shape become broader as the length of the free region w becomes 
large. It is worthwhile to mention that away from the resonance the value of Tr « is 
independent of the length of the well (see Fig. 12) and depends only on the barrier 
strength. The resonance being dependent on the incident energy of the electrons 
we follow a particular resonance peak (say, the first peak) in r^s for different E. In 
Fig. 13(6), the first peak values are plotted as a function of w for seven different 
incident energies in the sub-barrier tunneling regime. We see that as we increase 
the incident energy the resonances shift towards lower w value and also the peak- 
value decreases as expected. The above result clearly indicate that at and around 
the resonance the saturated value depends on the length of the free space, however, 
away from the resonance the saturated delay time is independent of the length of 
the free space ( space collapse or space destroyer as discussed earlier). 
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(a) 10 




(b) 45 



Fig. 13. (a) In absence of magnetic flux {i.e. <t> =0 ), for a ring with two identical barriers 
(Vi = V2 = 2, lb\ = lb2 = 5) separated by an intermediate well region of length w, the phase 
time Tr is plotted as a function oi w at E = 1.5. (&) With change in incident energy the position 
of resonance peaks shift. We follow the first resonance peak and plot resonant phase time r* as 
a function of w as we vary incident energy. Prom left to right the points correspond to energies 
E = 1.6, 1.5, 1.4, 1.3, 1.2, 1.1, 1.01 respectively. 



Instead of a real barrier we now consider a complex potential (or absorption) 
V = Vre — iVim'^^ Fig. 9(a). Here, Vre represents the real strength of the barrier 
whereas Vim is the imaginary strength. In presence of the imaginary part in the 
potential our ring system can absorb particles in it. Earlier Raciti et. al., have 
shown**^ that the phase time in transmission mode for an absorptive one-dimcnsional 
rectangular barrier does not show the saturation in opaque barrier limit, instead, 
it is linearly proportional to the length of the barrier. Thus they conclude that in 
presence of strong absorption the; tunneling time does not show the Hartman effect. 
We have separately stiidicd the phase time in both transmission and reflection 
modes of the one-dimensional complex barrier. We have found that in presence 
of non-negligible absorption in the large length limit of the barrier though the 
transmission phase time docs not saturate the reflection phase time shows Hartman 
eflfect. It would be worth to mention here that the saturated reflection phase time of 
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Fig. 14. In presence of a complex barrier covering the ring's circumference the phase time Tr in 
reflection mode is plotted as a function of ring's circumference. In (a) the three different curves 
are for three different strength of the real part of the potential while the imaginary part is fixed at 
Vim = 10- The solid, dashed and dotted curves arc for Vre = 2,4 and 10 respectively. In (6) the Tr 
is plotted for two different values of the imaginary part of the potential keeping the strength fixed 
at Vre = 2. The solid and dashed curves are for Vim = 5 and 10 respectively. Incident energy is 
set to be E = 1. 



ID complex barrier is same as that for an one-dimensional complex step potential. 
Now we study the phase time in reflection mode (Fig. 9(o)) in the presence of 
a complex potential. We vary the barrier length (ring circumference) for different 

values of Vre and Vim- From Fig. 14. wc sec that even in presence of complex 
potential Hartman effect survives for our ring system. In Fig. 14(a), keeping the 
imaginary part of the potential fixed we vary the strength of the barrier. We see 
that the saturated phase time as well as the nature of evolution of the r,. as a 
function of L depend on Vr,,- From Fig. 14(6), we see that the saturated reflection 
phase time changes with Vim as well. 

Finally we consider a quantum ring connected with two identical leads as shown 
in Fig. 15. Each of the two arms (say upper and lower arm) of the ring arc cov- 
ered entirely by a barrier. We study the transmission phase time in the absence 
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Fig. 15. Schematic diagram of a ring connected to two leads at the junctions Ji and J2 in the 
presence of an Aharanov-Bohm flux, </>. The upper arm of the ring is covered by a barrier of 
strength Vi and length Ibi. Similarly, the lower arm is covered by a barrier of strength V2 and 
length lb2 ■ The circumference of the ring is L = Ibi + i&2 ■ 

of magnetic fliox by varying different system parameters. Here also the phase time 
evolves with increasing length of the barriers and for larger length the phase time 
gets saturated. The saturation value does not depend on the ratio of the length of 
the two arms of the ring. However, this value depends on the incident energy of the 
particles. In Fig. 16 we have plotted the saturated transmission phase time Tts as a 
function of incident energy E/V for (f> = 0. The circumference of the ring is taken 
as i = 30 with equal upper and lower arm lengths. Plots with different arm length 
ratios {Ibi : ^62) with different (p in the asymptotic limit were found to overlap with 
the above curve in the entire energy regime. Analytically, in the large L (> 1/k) 
limit, the transmission phase time t/ becomes independent of L and the magnetic 
flux (in accordance with Hartman effect) and is given by 

2A;((2k2 - (fc2/2))2 + 4fc2K2)- ^ ' 



4. Conclusions 

We have verified the Hartman effect in quantum network consisting of an one di- 
mensional arm having two side branches. These side branches may or may not have 
barriers. In presence of a barrier, the 'phase time' for transmission through a side 
branch shows the 'Hartman effect'. Due to quantum nonlocality the 'phase time' 
and its saturated value at any side branch feel the presence of barriers in other 
branches. Thus one can tune the saturation vahie of 'phase time' and consequently 
the superluminal speed in one branch by changing barrier parameters in any other 
branch which is spatially separated from the former. Moreover, Hartman effect with 
negative saturation times (time advancement) has been observed for some parame- 
ter values. The phase times are also sensitive to the junction S'-matrix elements used 
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Fig. 16. In absence of magnetic flux (i.e. 0=0), for a ring shown in Fig. 15 the saturated trans- 
mission phase time rt s is plotted as a function of incident energy E/V. The length of the identical 
barriers (Vi = V2, Ibi = lb2) in the upper and lower arms of the ring are L/2 = 15. 

for a given problem. Depending on Wn there may be one or several bound states 
located between the barriers in different branched arms and as a consequence satu- 
rated delay time can be varied from the negative (ultra Hartman effect) to positive 
and vice-versa. We have verified this by looking at the transmission coefficient in 
the second arm ^2 which exhibits clear resonances as a function of W2 ■ Moreover the 
reported effects are amenable to experimental verifications in the electromagnetic 
wave-guide networks. 

We have further extended our studies on Hartman effect in quantum ring geom- 
etry in the presence of AB flux. We have studied the quantum reflection phase time 
for a ring connected with a single lead. Our study show that the phase time for a 
given incident energy becomes independent of the barrier thickness as well as the 
magnitude of the flux in the limit of opaque barrier. In the absence of AB-flux we 
have obtained analytical expressions for the saturated reflection and transmission 
phase times for the ring geometries. In addition, introducing a potential well be- 
tween two successive opaque barriers covering the entire circtimfercnce of the ring, 
we have found that the saturated reflection phase time becomes independent of the 
length of the well for energies away from the resonances. This implies, as if, the 
effective velocity of the electron within the well becomes infinite or equivalcntly 
the length of the well does not count (space collapse or space destroyer). Earlier 
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studies reveal that the Hartman effect disappears in presence of complex barrier 
in a transmission mode. We have shown that the Hartman effect survives in the 
reflection mode. 
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